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Abstract We present a linear perturbation analysis of 
the magnetorotational instability in the presence of the 
cosmic rays. Dynamical effects of the cosmic rays are 
considered by a fluid description and the diffusion of 
cosmic rays is only along the magnetic field lines. We 
show an enhancement in the growth rate of the unsta- 
ble mode because of the existence of cosmic rays. But 
as the diffusion of cosmic rays increases, we see that the 
growth rate decreases. Thus, cosmic rays have a desta- 
bilizing role in the magnetorotational instability of the 
accretion discs. 
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1 Introduction 

Understanding the true nature of accretion processes 
in astrophysics has always been an attractive research 
topic over the last three decades. Accretion discs are 
observed in many astrophysical systems from new born 
stars to compact objects or even very large discs at the 
center of the galaxies. In spite of the diversity of the 
accreting systems, existence of a possible mechanism 
of the angular momentum transport is a common fea- 
ture in all these accretion systems. Extensive efforts 
to understand mechanisms of the angular momentum 
transport in the accretion discs have lead to a better 
understanding of such systems, though there are many 
theoretical and observational uncertainties. 

It has been proposed that the magnetorotational 
instability (MRI) is the main driving mechanism of 
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turbu lence in the accretion discs (jBalbus and Hawlev 
199lh . Extensive subsequent works have clarified and 
extended our understanding of the role of MRI in var- 
ious astrophysical systems fr om protoplanetary discs 
(e.g., Sano and Mivama 19991) to the protostellar discs 
or even quasar discs. Over recent years a multi layer 
model for the protoplanetary discs is proposed in which 
the surface layers are magnetically active due to the ion- 
ization of the CRs, while the central layers are magnet- 
ically inactive because of inabili ty of CRs to p enetrate 



down to the central parts (e.g.. iGammid 120011 ). Thus, 
MRI can act as a driving mechanism of the turbulence 
and the accretion at the surface layer in a protoplane- 
tary disc. However, possible dynamical effects of CRs 
on MRI have not been studied to our knowledge. 

Cosmic Rays are very energetic particles but their 
energy density is in equipartition wit h energy dens i- 
ties of thermal gas and turbulence (e.g., FerrierelboOl ). 
CRs can act as a source of heating and increase the 
level of the ionization in the interst ellar medium (e.g 



Field. Goldsmith, and Habing||l969l) . An enhanced flux 
of CRs has important con sequences for star formation 
near to the Galactic center ( Yusef-Zadeh. Wardle. and Rov^ 



20071) 



However, interaction of CRs with a plasma is not 
restricted just to a possible enhancement of the level 
of ionization as have been studied extensively over re- 
cent decades. For example, dynamical effect of CRs 
has a vital role in analysis of Parker instability for the 
structu r e formation in the Galaxy a.t large scale (e.g., 

Parker 1966 : MouschoviasI 1974 : ^Hanasz and Lesch 

200C ; Kuwabara. Nakamura. and Kcii2004i : iKuwabara and Ko^ 



20061) . It is also found that C Rs have a stabilizing ef- 
fect on the thermal instability (Kuznetsov and PtuskinI 
1983t IWagner et aZ.ll2005l : IShadmehrill2009l ). 

The problem of the diffusion of CRs and its role in 
MRI has not been studied in detail. Considering CRs as 
a separate fluid and their diffusion along the magnetic 
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field lines, we study MRI in the presence of CRs via a 
linear perturbation analysis. Our basic equations and 
the assumptions are presented in the next section. Final 
dispersion relation will be analyzed in sections 3 and 4. 



2 General Formulation 

In this study, CRs are protons, electrons and nuclei. 
But we neglect the electrons because of their little con- 
tribution to the total pressure. There are three dif- 
ferent approaches to study the dynamics of CRs. In 
the particle-particle approach, the plasma and CRs are 
considered as particles that may interact with each 
other via complicated processes. In a simpler approach, 
known as fluid-particle, the plasma is treated as a fluid, 
though CRs are still described as particles. The sim- 
plest approach is the fluid-fluid approach in which CRs 
and the thermal gas are described by different interact- 
ing fluids. The hydrodynamic approach can not provide 
the spectrum of CRs, however, it is a good approxima- 
tion for analyzing dyna i nics of a plas ma with CRs (e.g., 
Drurv and Voelklll98ll IPrurvlllQSsI) . Thus, we adopt 
the hydrodynamic approach to study the effects of CRs 
on the unstable modes of MRI. 

We also follow the same steps as in, except that CRs 
arc considered as a fluid, and diffusion is considered 
only along magnetic field lines. For simplicity, we ne- 
glect the ionization and the heating by the CRs, since 
their effects in the absence of the dynamical role of CRs 
are well understood. 

The basic equations are 



dp „ 



= ^(V X B) X B- V(p+Pcr), 

dB 



dt 



V X (v X B), 



(1) 
(2) 
(3) 



field, respectively. Also, the net cooling function is de- 
noted by A in the energy equation ([5]). However, we 
will neglect the net cooling in this paper, i.e. A = 0. 
Moreover, we have 



-K||b(b • Vpcr), 



(7) 



is the diffusive flux of cosmic-ray energy and Ky is dif- 
fusion coefficient along magnetic fleld lines. All the 
variables have their usual meaning. Also, b is a unit 
vector along the magnetic field lines, i.e. b = B/i?. 
The adiabatic indices of the thermal gas and cosmic 
rays are denoted by 7g and 7cr, respectively. Finally, 
we can write equation of state as 



P 



-pT, 
P 



(8) 



where R is the gas constant and fi represents the molec- 
ular weight. 



3 Linear Perturbations 

Assuming that the system is axisymmetric, we write 
basic equations in the cylindrical coordinates (i?, <&, z). 
The equilibrium magnetic field is assumed to be con- 
stant in space and only the toroidal and the vertical 
components of the magnetic field are considered, i.e. 
Bq = -Bz^z + i?$e$. The equilibrium disc is rotating 
with Keplerian angular velocity, i.e. il ~ ^ GMj B? 
where M is the mass of the central object. Moreover, 
the initial density and the gas and the cosmic ray pres- 
sures are considered to be constant. Now, we can per- 
turb the equations of the form bX ~ SXQCXTpi{k^R + 
kzZ — ujt). Thus, the linearized equations become 



S p 

'iuj h ifcR(5wR -I- ikzSvz = 0, 

P 

-iujSvYi — 2r2(5w$ + —^Sp + —^Spc 
P P 



(9) 



V • B = 0, 



de 



P-^ = -PV • V A, 
and 



1 dpc 



7cr Per dp 



7cr - 1 dt 7cr - 1 p dt 



V - r = 0, 



(4) 



(5) 



(6) 



where d/dt ~ d/dt+v-W is the Lagrangian time deriva- 
tive. Here, p, p, pcr, v and B are the density, gas pres- 
sure, cosmic ray pressure, velocity and the magnetic 



-^BJBn + ^{B^SB^ + BJB,) ^0, (10) 



-iujSv^ + -^tttSvu - —^BJBis. = 0, (11) 
2,\l Airp 



.,^5^^ + '^5p + 'J^5p^^ + '^BeSBe = 0, (12) 
p p Anp 



-iujSBji ~ ikzBzdvB. = 0, (13) 
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-ikzB^Svz — 0, 



dlnR 



(5 -Br + iki{Bq,Svji 



(14) 



Z^CoQ!^(7g + 7cr'P)] cose 

Pe = 2 7g Co ■00 9 "^^^^ 2'(7cr - 1) cos^ 6 



(25) 
(26) 



-iujSBz + ik^^BzSvB. = 0, 



5p Sp 

-iuJ h iuj^e: 

P P 



7c 



1 



K||(b.k)2 



Spc 



7c 



Pc 



7c 



1 p 



where is the square of the epicycHc frequency, 



(15) 



(16) 



■iiu6p = 0, (17) 



(18) 



Introducing a new dimensionless variable as X 
uj/{iD,), our final dispersion equation becomes 



P7X^ + PqX^ 



PaX^ + P:<X^ 



P2X^ + PiX + Po = 0, 



(19) 



where the coefficients arc complicated function of the 
input parameters. We denote the angle between the 
vector k and the component by Q. Also, a nondimcn- 
sional wavenunibcr is defined as Z = kzVAz/^- Thus, 
the coefficients become 



Po = 2 7g A V"? a'f^'co'Z'^ (7er - 1), 



(20) 



P7 = 2 7g a^n^ cos'^e 



(27) 



where A = 2A cos^ Q + Z^. We verified that the above 
dispersion equation reduces to the equation for a case 
without CRs (e.g., see Eq. (17) of Sano and Mivamal 
( 1999! )). In the above relations, we have A = — (i(lnJ7)/d(lnP),| 
q = [b.(k/fc)]^. Also, the parameter (p represents the 
ratio of the pressure of CRs to the gas pressure, i.e. 
(f = Pci/p- The initial direction of the magnetic field 
is denoted by a = B^/B^. The plasma beta param- 
eter is defined by the poloidal field, /3z = 2c'^/-fgvj^^ 
where Cs is the sound speed and vaz is the Alfven 
velocity, i.e. uaz ~ Pz/V^T^p. Thus, we obtain 
Cs = cqwAz, where cq = ^ (7g/2)/3z. Also, we have 
ip = K||fc/cs = ipoZ/ cosQ, where the nondimensional 
diffusion coefficient ipo is 



00 = Co 



TD 



(28) 



where the diffusion time scale td, the dynamical time 
scale Td and the sound crossing time scale Tc are defined 
as 



P^ 



1 



P 



,Td = o,Tc = — . 

K\\ il Cs 



(29) 



Pi = 2 A a^n^ co^ Z"*(7g + 7cr 0)cose, 

P2 = 2[(1 + cla^){K'^/n^) cos^ e + a^A + (1+ 



2cla^)Z^]{^„-l)n^Z^jgCotPq 

P3 = -Z^n^ cos^ e{7gCoV'qZa(7cr - 1)(4 + 2A 

-kV^^^) sine - 2cos2 e[cga2 (7g + 7cr<y3)KV^^^ 
+7g(2Aa2 + nyn^)] - 2Z^[2cla^{jg + 



(21) 



(22) 



(23) 



+7g(l + a')]} 
P4 = -/f,Zn^[a{n'^/n'^ - 2A - 4)Zcos2 esine + 2coX 
0'oga^(7cr - l)iZ^ - K^/n^) cos^ e + 2co0og(7cr - l)x 

{l + a^'+cla^)Z^] (24) 
P5 = 2n^['jga^{Z'^ + K^/Vl?) cos^ e + 7gZ2(l + a^)+ 



4 analysis 

Equation describes magnetorotational instability 
with CRs. Considering complexity of the coefficients of 
our dispersion equation, it is very unlikely to obtain the 
roots in analytical closed forms. However, we can solve 
the equation numerically for a disc with the Keplcrian 
angular velocity. Since our goal is to analyze the un- 
stable perturbations, we restrict our study to the roots 
with positive imaginary part, i.e. Im(aj) > 0. If we ne- 
glect terms corresponding to the CRs, i.e. = ■00 = 0, 
our algebraic dispersion relation (fl^l) reduces to the 
classical dispersion equation (e.g., ISano and Mivama 
I999I) . Among our input parameters, the effects of CRs 
are described via two input parameters and "00. One 
can vary these parameters to study how the unstable 
modes are modified due to the existence of CRs. 

Typical influence of CR on the magnetorotational 
instability are shown in all subsequent plots. Fig- 
ure [1] shows nondimesnional growth rate of the unsta- 
ble mode, ljj/Q., versus nondimensional wavenumber Z. 
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Fig. 1 Imaginary part of the growtli rate as a function of tlie vertical wave number 
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Each curve is labeled by its corresponding parameter. 
We found the maximum growth rate occurs for the per- 
tur bations with fcR, = (see also, Sano and Mivamal 
[l999il . The top left-hand plot of Figure [T] shows behav- 
ior of the unstable perturbation when the ratio of the 
CRs pressure to the gas pressure changes from zero 
to one with ■^o = 0, = 0, /3z = 20 and a = 0.1. Here, 
diffusion of CRs is neglected. Obviously, the case with 
iy9 = corresponds to the unstable mode without CRs. 
As the ratio Lp increases the system becomes more un- 
stable because of the enhancement of the growth rate. 
In other words, existence of the CRs destabilizes the 
disc. 

The top right-hand plot of Figure [1] shows the typ- 
ical dependence of the growth rate on the parameter 
Pz which is ratio of the gas pressure to the magnetic 
pressure. The input parameters are V'o = 0, = 0, 

= 0.6 and a = 0.1 and 13 z varies from low value 0.1 
to high value 100. Corresponding to each case repre- 
sented by the solid lines, there is a dashed curve which 
is for the same case but without CRs. The bottom left- 
hand plot of Figure [U shows the case with CRs (solid 
line) and without CRs (dashed line) for different values 
of a which is the ratio of the z component of magnetic 
filed to the ip component of the magnetic filed. The 
parameter a changes from low value 0.01 to high value 
0.5 and the other input parameters are ^/iq = 0, = 0, 
Lp = 0.6 and f3z — 20. In all previous plots, diffusion of 
CRs along the magnetic field lines is neglected. Now, 
the bottom right-hand plot of Figure [1] shows how the 
growth rate of the unstable perturbation is modified 
when the diffusion of CRs is not negligible. Here, all the 
input parameters arc fixed asa = 0.1, = 0, Lp — Q.Q 
and (3z = 20 but the dimensionless diffusion parame- 
ter 'i/'o changes from to 1000. Diffusion of CRs has a 
stabilizing effect according to this plot. However, when 
diffusion is high we can hardly recognize any changes 
and the plots overlap. 



along the magnetic field lines increases, CRs pressure 
decreases and the system tends to a case without CRs. 
Thus, it would be interesting to study MRI with CRs 
at the nonlinear regime by doing numerical simulations. 
There are regions with high fiux of CRs such as near to 
the Galactic center or ultra luminous infrared galaxies 
( PapadopouloslliolO : Papadopoulos et aZ.|[201ol) . Under 
such circumstances, we think, evolution of the accretion 
discs located at such regions are significantly affected 
by the dynamical effects of CRs. 
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5 conclusion 

Our simple approach shows that CRs destabilizes MRI 
unstable modes. It implies that the generated turbu- 
lence because of the MRI would be amplified in the 
presence of the CRs. However, it seems that growth 
rate only slightly modifies in the presence of CRs in 
the linear regime. In our model, gas and CRs coupling 
is via cosmic ray pressure term in the equation of mo- 
tion. In other words, CRs provides an extra pressure 
and we know the pressure plays minor roles in MRI so 
long as the magnetic pressure is smaller than the gas 
and the CRs pressures. But when the diffusion of CRs 
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